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1 Introduction 

The finite volume element (FVE) method has been widely used in numerically solv¬ 
ing partial deferential equations. The main feature of FVE method is that it inherits 
some physical conservation laws of original problems locally, which are very desir¬ 
able in practical applications. During the last decades, many research works have 
been presented for FVE methods solving various partial differential equations, see 
lfni2l[^l4ll5ll6ir71l^l9i rT()ll 1 211 1 3ll 1 411 1 SIITtII 1 8II20I and the references cited therein. 

Superconvergence of numerical solutions has been an active research area for 
finite element method (FEM) since its practical importance in enhancing the accuracy 
of finite element approximation. But the study of superconvergence properties of FVE 
methods is far behind that of FEMs. For elliptic problems in two-dimensional domain, 
the early superconvergence results of linear and bilinear FVE solutions are of this 
form imiTOl 

<Ch^\\uh^^, ( 1 . 1 ) 

zes 

where S is the set of optimal stress points of interpolation function on partition 7},, 
N = 0{h^^) is the total number of points in S, and V denotes the average gradient 
on elements containing point z. Eater, Lv and Li in ifT^ extended result (11.11) to the 
isoparametric bilinear FVE on quadrilateral meshes under the fi^-uniform mesh con¬ 
dition. Recently, Zhang and Zou in Il20l also derived some superconvergence results 
for the bi-complete k-order FVE on rectangular meshes, and in the case of bilinear 
FVE {k= 1), their result is 

\Wiu-Uh)iG)\<Ch^\lnh\^\\u\U,^, ( 1 . 2 ) 

where G is the Gauss point of element (the midpoint of element). Moreover, by using 
the postprocessing technique, Chou and Ye El obtain the superconvergence estimate: 

||Vm-V gM/,)|| <C/!2||m||3, 


(1.3) 
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where Qu;, is the postprocessed linear FVE solution obtained by the L 2 -projection 
method; Zhang and Sheng IfTSl further derive the superconvergence estimate: 

||VM-f;VM;,)|| <C/i2||m||3, (1.4) 


where is the reconstructed gradient of bilinear FVE solution obtained by using 
the patch interpolation recovery method. 

In this paper, we consider the bilinear FVE method to solve the following problem 


—div{AVu)-\-cu = f, in Q,, 
M = 0, on dQ, 


(1.5) 


where C is a rectangular domain with boundary dQ, coefficient matrix A = 
{aij) 2 x 2 - Our main goal is to give some piecewise-point gradient superconvergence 
for the bilinear FVE approximation to problem (ll.51 l. To the authors’ best knowledge, 
in existing literatures, only the midpoints of elements are proved to be the supercon¬ 
vergence point of gradient approximation IMI (also see (fO l. we here will prove 
that except the midpoints of elements, all interior mesh points and midpoints of inte¬ 
rior edges are also the superconvergence points. Generally speaking, the analysis of 
bilinear FVE on rectangular meshes is more difficult then that of linear FVE on trian¬ 
gle meshes, the reason is that Vm/, is not constant in the former case. By calculating 
exactly some integrals on element and its boundary, we first establish the superclose 
weak estimate for the bilinear form of the FVE method. 


<C/!^||M|| 3 ,p||v||r 5 , VvG t//,, 2<p < oo, l/p+l/q= 1, (1.6) 


where fl/j is the bilinear interpolation operator, 17^* is the projection operator from 
the trial function space t//, to the test function space. It is well known that such weak 
estimate plays an important role in the superconvergence analysis of FEMs 111 1II19I 
ED. By means of this weak estimate and some integral estimates on element, we 
further derive the following superconvergence result for the gradient approximation 

max |(Vm — Vm/j)(P)| < C/z^j ln/!|||M|| 3 ^cx,, 


(1.7) 
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where S is the optimal stress point set of bilinear interpolation function which is com¬ 
posed of the interior mesh points, the midpoints of interior edges and the midpoints 
of elements. 

This paper is organized as follows. In Section 2, we introduce the FVE method 
and give some lemmas. In Section 3, we establish the superconvergence weak esti¬ 
mate. Section 4 is devoted to the piecewise-point gradient supconvergence estimate 
on the set S. Finally, in Section 5, numerical experiments are provided to illustrate 
our theoretical analysis. 

Throughout this paper, we adopt the notations W'”’P{D) to indicate the usual 
Sobolev spaces on domain D C Q equipped with the norm || • \\m,p,D and semi-norm 



we omit the index D. We will use letter C to represent a generic positive constant, 
independent of the mesh size h. 

2 Finite volume element method 

Consider problem (11.5b . As usual, we assume that there exist positive constants Ci 
and C 2 such that 


GR^, (x,y)GQ, 


( 2 . 1 ) 


We further assume that A G c G L„{Q,) and c > 0. 

Let 7), = U{^} be a rectangular partition of domain Q so that Q. = U/rer;, { ^ 
where h = maxhK, hx is the diameter of element K. We assume that partition T/, is 
regular, that is, there exists a positive constant 7 > 0 such that 


hx/PK < £ 7/2) 


( 2 . 2 ) 


where px denotes the diameter of the biggest ball included in K. 

Associated with partition 7),, we construct the central dual partition 7)* by con¬ 
necting the center of each element to the midpoints of edges by straight lines. Thus, 
for each nodal point P in T/,, there exists a rectangle Kp surrounding P, Kp G Tj* is 
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called the dual element or the control volume at point P, see Fig. 1. 





p 










FIG. 1. The dual element Kp sun'ounding point P 

To partition 7/, and 7}*, we associate the following trial function space 77/, and test 
function space Vk, respectively, 

Uh = {uh€C%n)-.Uh\K€Qx{K),^K€Th, M/,Ui2=0}, 

Vh = {v/, G L2{Q.) ; Vh\K*p = constant yp G Nh, = 0, VP G dPl}, 

where Q\(K) is the set of all bilinear polynomials on K and Nh is the set of all nodal 
points of Th- 

The conventional weak form for problem (11.51) is that find m G 77q (X2) such that 


a{u,v) = (/,v), VvG77o'(f2), 


(2.3) 


where 


(2.4) 


a{u,v) = / AVm-Vv + cmv, (/,v)= / fv. 

Jn Jn 

Weak form (12.31 ) usually is adopted for the finite element method. But for the FVE 
method, we need a new weak form. Let u be the solution of problem (11.51 ). then by 
using Green’s formula, we have 


■ / n-{A'Vu)vds+ / CUV = / fv, 
JdK*„ JKi JKt 


K*p G Tf, V G Vh, 


(2.5) 


where n is the outward unit normal vector on the boundary concerned. Motivated by 
this weak formula, we introduce the following bilinear form 

«/,(«, v)= F ( — / n-{A'Vu)vds+ / cuv\ u G {£2), v GVh, (2.6) 

K*eT* ^ •'^p ^ 

and define the FVE approximation of problem (11.51) by finding m /, G Uh such that 


ahiuhNh) = {f,Vh), Vv/, G Vh- 


( 2 . 7 ) 
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Let : Uh Vh be the interpolation operator defined by 

n*,vf,= Y, vhiP)Xp,^vh ^Uh, 

PeNi, 

where Xp is the characteristic function of the dual element Kp. Since is a one to 
one mapping from t//, onto V/,, then the equivalent form of problem (I2.71 i is that find 
M/i G Uij such that 

ahiuh,n^,Vh) = Vv/, G Uh- (2.8) 

This is the FVE scheme to be used in our analysis. From (I2.5l l we know that scheme 
( 12 . 8 b is consistent and the following error equation holds. 

ah{u-Uh,n^Vh)=Q,yvhGUh. (2.9) 

Let fl/jM be the usual bilinear interpolation function of continuous function u. 
In our analysis, the following approximation property, trace inequality and inverse 
inequality will be used frequently. For 1 < p,q <°o,we have 

\\u-nhu\\m^p^K<Ch^f"'\\u\\2,p,K,0<m<2, (2.10) 

\W\\a,p,dK < Chpr'’ (||M||o.p,/t + %||VM||o,p,^:), u G (2.11) 

\\uh\\m,p,K<Chl^~^h'^’”\\uh\\i^^^^ UhGQx{K),Q<l<m<\. (2.12) 

Furthermore, the following two lemmas hold II10II18I . 

Lemma 2.1 Let K G T),, t C dK be an edge ofK. Then, for Vh G Up, I < q < °°, we 
have 

[ (y, -= 0, /(v, - n^Vh)ds = 0, (2.13) 

Jk Jr 

\\n-n*pVh\\o,,,K < ChK\\vh\\i^,^K, (2.14) 

H-n*pvp\\o^p^3K<ch^^h\nhq,K- (2.15) 

Lemma 2.2 For h small, we have 

ah{vh,nhVh) > C\\vh\\], Vv/, G Up. 

Lemma 2.2 implies that the solution m/, of problem (12.8b uniquely exists. 
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3 Interpolation weak estimate 

It is well known that the interpolation weak estimate plays an important role in the su¬ 
perconvergence analysis 01 llll9ll2T1l . To FEM defined on regular rectangular meshes, 
the following interpolation weak estimate has been established 01911211 

\a{u- nhU,v) \ < C/!^||m||3,p||v||u5, VvG 1//„ 2 <p < oo, 1/p+ 1/^ = 1. (3.1) 

In this section, we will give a similar weak estimate for the EVE method, which is 
very useful in our superconvergence analysis. 

In order to utilize the known result (13.10 . we need to give the difference between 
the bilinear form of EVE and that of EEM. 

Lemma 3.1 For any w £Uit+H^{i2),Vh £ Uh, we have 

ah{w,n^Vh)-a{w,Vh) = T / n-{AVw){n^Vh-Vh)ds 

K€Ti, 

+ {-div{AVw) +cw,nfjVi,-Vh)K, (3.2) 

Ken 

where Uh + FP' (12) = {w = «/, + v : «/, G Uh,v € [Q .)} is the the algebraic sum 

space. 

Proof By Green’s formula, we have 

/ A'^w ■'Vvh = — / c//v(AVw)v/, + / n - {AVw)vhds, 

Jk Jk JdK 

and (see Eig.l) 

Y [ div{AVw)nhVh = Y H i div{AVw)nhVh 

Ken-'^ Ken K*eT* 

= Y [ n-{A'Vw)nhVhds+ Y [ ^n- {A'Vw)nhVhds. 

Ken KpeT^ ''^^*p 

Substituting this two identities into the definitions of a{w,Vh) and ah{w,nf^Vh) (see 
(1240 and (l231 l). the proof is completed. □ 

Let/T = \I\PiP 2 Pj,Pa be a rectangular element, hx = X 2 — xi, hy =y2 — yi, see Eig.2. 
The following lemma gives some exact calculations for integrals on dK. 
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P4(Xry2) 


M 


Pi(Xi,y, 


P3(X2,y2) 


P2(X2.y,) 


FIG.2. Rectangular element K — □PiR2^3^4 
Lemma 3.2 For w, v € Uh, we have 

hi 


[ in^,v-v)wxdy= ^Vy{xi,y)wjcy, f {n^,v-v)wydy = 0, (3.3) 

JP 1 P 4 24 Jpj/>4 


[ i^hV- v)wxdy = y^Vy{x 2 ,y)wjcy, f (!//>- v)wydy = 0 , (3.^ 

JP2P3 24 JP2P3 


f h^ f 

/ inhV-v)wydx= 2 :^Vxix,yi)wjcy, / [U^v-v)wxdx = 0, (3. 

JP1P2 ' JPiPi 


4) 


5) 


6 ) 


f f 

/ inhV-v)wydx= y:^Vx{x,y 2 )wxy, / {n^v-v)wxdx = 0. (3. 

JPyPy 

Proof We only prove (I3.31 I. the other are similar. Noting that Wx and v are linear on 
line segment P 1 P 4 (M is its middle point), we have from the definition of TI^ v that 


/ n^vwxdy = / v{Pi)wxdy+ v{Pf)Wxdy 
JPiM JMPy 


= ^v{Pi){wxiP\) +Wx{M)) + ^v{Pa){wx{M) +WxiPA)) 
= ^{v{P,)wx{Pi)+2v{M)wx{M) + v{P^)wx{Pa))- 


On the other hand, using the composite trapezoidal rule: 
f(v')iiv = - ( f(o\ A-7.f( — ^ 


rb 

/ f{y)dy - 

J a 




we obtain by taking / = vwx that 
h 

vwxdy = -^{v{P\)wx(P\) +2v{M)wx(M) PvfPfjWxiPA}) - ^Vy(vi,y)wxy. 

The difference of above two equalities implies the first equality in (I3.31 I. The second 
equality comes from (I2.13I) . noting that Wy is constant on P 1 P 4 . □ 
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Let denote the piecewise constant approximation of function w on 7/,, 

w‘'\k = ^ [ w, K GTh; \\w-w‘'\\o,p,K<ChK\\w\\i^p.K, (3.7) 

Theorem 3.1 Let 7), be a regular rectangular mesh and u € Then we have 

|a/,(M-7T/,M,77;»| <C/!^||m||3,pIM|i,^, WvGUh, 2 < p < oo, l/p+ l/^= 1. (3.8) 

Proof Let Am = (dij) 2 x 2 be the value of A at the midpoint M of edge T C dK, see 
Fig.2. From Lemma 3.1, we have 

ap{u — TlfiUjll^v) —a{u — TlhU^v) 

= y f n-{A-AM)'^iu-n,,u){n^v-v)ds 
Ken 

+ y n-{AM'^{u — nhu)){nlv — v)ds 
Ken 

+ ^ {-div(A'V{u-nhu)),n^v-v)K 
Ken 

+ E {c{u - nhu),niv- v)k 
Ken 

= £i +£2 + £3+£4. (3.9) 

Using (I2.10I) - (I2.11I) and (I2.14I) - (I2.15I) . we obtain 

£1 < C ^ hK\A\v^\\W{u-nhu)\\o,p^dK\\v-n*^^v\\Q^q^3K<Ch^\\u\\2,p\Muq- 
Ken 

£4 < c ^ ||M-7T/,M||i.p,;,||v-77>||o,,,ir<Ch2||M||2,p||v||i.,. 

Ken 

For £3, setA = (31,82). Since 

div{A'Ww) = (^7^val,^7^va2) • Vw + ai-Vwj: + a 2 -Vwj,, (3.10) 

then we have 

£3 = ^ {—{divai^divaz) ■'^{u — nhu),n^v — v)K 

Ken 

- (ai-V(M-77pM);t + a2-V(M-77/,M)j.,77,*v-v)/f = £3 i+£32. 

Ken 
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Obviously, £31 < C/z2||m||2,p||v|| 1,5. Now, let «,• = u{Pi), ^ = {x-x\)/h^, 77 = (y- 
y\)/hy, see Fig. 2 . Since 

UfiU = Ml + (m 2 — U\)^ + (m 4 — Ml)?] + (m 3 + Ml — M 2 — M 4 )^ 77 , (.x,y) S 


hence 


{ni,u)xy = (m3 + Ml - M2 - M4)/!^ ' 

1 /■j :2 m 

=-p\ Ux{x,y2)dx- Ux{x,yi)dx) 

A Jxy Jxi 

1 m ryi l r 

= — / / Uxydxdy = — Ujcydxdy = {u^f. 

Jx Jx, Jvi £ 


■'Jti ■'yi 

Therefore, we have from (12.131 1 that 


(3.11) 


£32 = - ^ (aiiM;cA: + a2i(Mxy-(n/,M)j^)+ai2(Mxy-(I//,M)^)+a22Mrr,F/p*v-v)A: 
Ken 

= - X (ailM« + a2l(Mxy-«Jy) + ai2(Mj3--Mjy)+a22M>")-,npV-v)A: 

Ken 

= - Yj {a\\Uxx + a22Uyy-{aiiUxx + a22Uyyy,nf^v-v)K 

Ken 

- Y i‘^2l{Uxy-u‘xy)+an{Uxy-u‘xy),n^V-v)K- 

Ken 

Using the approximation property, we obtain £32 < C/i^||m|| 3 ,pIM|p^ so that 


£3 = £31 +£32 < C/!^||m|| 3 _p||v|| 


Now we need to estimate 

E2=Y Y f n-{AM^{u-nhu)){n;v-v)ds, (3.12) 

Ken tcdK\dn 

noting that {FI^v — v)|5^ =0. Let T be an interior edge, that is, a common edge of 
two adjacent element K and K'. Since nltn^A" = ~”lTn5A'' ^nd Am^ u{nf^v — v) is 
continuous across edge T, then we have 

E2=Y Y (-n-{AMVnhu){nlv-v)ds= Y Y (3.13) 

Ken xcdK\dn Ken TedK\dn 
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Without loss of generality, let K = = ixi,X 2 ) x (yi,y2) be a rectangular 

element and the common edge T = dK f] dK' = P 1 P 4 , see Fig.2. From Lemma 3.2 we 
have 

r 

F{z)= / {an{nhu)x + ai2{nhu)^.){n^v-v)dy = an::^Vy{xuy){nhu)jcy. (3.14) 

Then, from (13.141 1. (I3.4l i and (13.1 11 1, and noting that n|Tn3A: = 
is continuous acrose T = dK fl dK' = P 1 P 4 , we obtain 

F{^C^dK)+F{^C^dK')=au ^v^(xi ,y){{u^y{K) - {u^y{K')). (3.15) 

Let £ C®(X2), ihw\K & Qi (K), be an approximation of function w, which satisfies 

\\w-ihw\\o.p < Ch\\M\ i,p> 2 < p < 00. 

For example, for p — 2we may take ipw — Ppw, Pp is the L2-pi'ojection operator; For 
p > 2, we may take //,w = Tlpw. Then, from (13.15b and the inverse inequality, we have 

F{TC^dK)+F{xC^dK') 

_ h^, 

— ^ii'^^yixi:y)yuxy) {K) iiiUxy{M) iijUxy{M') {i^xy) {F y 

_2 _2 _2 
"^Chyh^ ||t'|| ihUxy\\0,p,K F h^, ||(M4;y) Ui^xy\\Q.p,Ky 

< C/jvllvjl (II (uxy) — Mjiyllo.p.ArjAT' + 11^431 — i/iMj 3 i||o,p,ArjAr') ^ ( 3 . 16 ) 

where we have used the fact that hy = \dKy\dK'\ < min{/r/i:,/i^/}. Combining (13.13b 
with (13.16b . we can conclude that 

£2 < Ch\\v\\l^q{\\Ujcy - {u^yh^p + IIma,, - 4Mxy||o,p) < C/i^I|m|| 3 ,pIM| l,q. (3.17) 

Substituting estimates E\ ^ £4 into (13.9b . we complete the proof by using (13.1b . □ 

From Theorem 3.1, we immediately obtain the following superclose result. 
Theorem 3.2 Let u and u/, be the solutions of problem (O and ( I2.1SI) . respectively, 
u £ H^{£2). Then we have 


WhU-HhWi < Ch^\\u\\2. 
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Proof From Lemma 2.2, error equation (I2.91 l and weak estimate (I3.81 l. we have 

< ah{uh - nhU,n*^[uh - Tlhu)) 

= «/,(« — — riltu)) < C/!^||M||3||M/i — i7/,M|| 1. 

This gives the conclusion. □ 

By using the embedding theory and weak embedding inequality of finite element 
space ED, a direct result from Theorem 3.2 is the following optimal Lp-error esti¬ 
mates, 

\\u-Uh\\o,p < ||m-I//,m||o^p+ ||fT/,M-M/,||o.p < Cp/z^(||M||2,p+ Hulls), I <p<°°, 
||m-m/2||o,oo < C(/!^||m||2.oo+ |/n/r|^||fT/,M-M/,111) < C/r^(||M||2.«,+ |ln/z|2||M||3). 

4 Piecewise-point superconvergence for gradient approximation 

Definition 4.1 Let point P G Q, we call P an optimal stress point for interpolation 
operator LIi,, if 

|(VM-VfT/,M)(P)| <C/z2||M||3.<»,i?, (4.1) 

where E is the union of elements containing point P and VwfP) denotes the average 
value of gradient Vw{P) in elements containing point P. 

It is well known that for the bilinear interpolation TlhU on rectangular meshes, 
the optimal stress point set S is composed of all interior mesh points, midpoints of 
interior edges and the midpoints of elements 016II21II . In this section, we will derive 
the VT*’“-superconvergence result for the bilinear finite volume solution m/, on set S. 

We first give some lemmas. Let rectangular element K = □PiPsPsA- For wp G Up, 
set Wi = w{Pi),i = 1, • • • ,4. Introduce the discrete //^norm by 

\wh\i.h = { Y, ^/2 G 14, (4.2) 

K€T^ 

where 

k/ill./i.A: = (>^2 -Wi)^ + (wj, - wff- + {W3 - wa)^ + (W4 -w\f-. ( 4 . 3 ) 
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Lemma 4.1 The discrete norm Iw/ili,/, is equivalent to the norm |w/,|i on space Uh, 
and 

^I'^/iIi,/!,a'< K/ilix^ VITS 7},. (4.4) 

Proof Let K = UPiP 2 Pj,Pa = [xi,xi + /tj x [yi,yi + hy] ,^ = {x-xi)/hx, r] = {y- 
y 1 ) /hy. For Wh & Q\{K),wt may write 

V^h{x,y) = Wi + (W2 — wi)^ + (W4 — Wi)77 + (w3 + Wi — W2 — W4)^r]. (4.5) 


Set 


1^21 = W2 — Wi, W41 = W4 — Wi, W32 = W3 — W2, W34 = W3 — W4, (4.6) 

W1234 = W3 + Wi — W2 — W4 = W34 — W21 = W32 — W41 . (4.7) 


Then, we have from (|43])-(l4j]l that 


\'>Vh\i.K = jj'^'^h\^dxdy = {{d^Wh/hxf+ {djjWh/hy)^)h^hyd^dr^ 

I 2 rl r\ ^ rl rl 

=-r / {w2\+wx2iAr]fd^dr] + -^ i / (w4i + wi234^)^<i^<7i? 
hx Jo Jo hy Jo Jo 


hy / 9 ^ 2 \ / 2 ^ 2 \ 

~ — (W 21 + W 21 W 1234 + ^^ 1234 ) ^41^1234 + 0 ^ 1234 ) 

hx 3 hy 3 

1 hy / 2 2 \ ^ / 2 2 \ 

~ O T^('^21 + VV 21 W 34 + W 34 ) + - —(W 32 + W 32 W 41 + W 42 ). 

3 hx 3 hy 

Hence, it follows from the Cauchy inequality and the regularity condition (12.21) that 


^(^21 +W34 + W32 + W41) < |w/,|i < 2 (w2i +W34 + W32 + W41). 


Combining this with (14.31 ). the proof is completed. □ 

Let = (aL) 2 x 2 be the piecewise constant approximation of matrix function A. 
Lemma 4.2 The following identity holds for w, v C Uh, 


f n-{A‘^Vw){nlv — v)d. 
JdK 


hlhx hyhl 

-V}as= -r^aiiVxyWxy+^^a 22 VxyWxy (4.8) 


24 
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Proof From Lemma 3.2 we obtain 

[ n ■ {A‘^'Vw)(n^v — v)ds 
JdK 

= / in^v-v){a‘liWx + a‘l2Wy)dy- / {n^v-v){a‘iyWx + a‘l2Wy)dy 
JP 2 P 3 JP 1 P 4 

+ / {n^jV — v){a2iWx + a22^y)dx— / {n^v — v){a2iWx + a22^y)‘^x 
JP 4 P 3 ' JPt Pi 


hi 


hi 


= ^a‘iivy{x2,y)wxy - ^a'iiVy{xuy)wxy 

h^ h^ 

+ :^a^22Vxix,y2)Wxy - ^a^22Mx,yi)Wxy 


24 

hl.hx 
y ^ 


hyhl 


24 ^ll^xyWxy+ 24 022VxyWxy- 


□ 


Now, we need to introduce the regularized Green function 019II21I . For any given 
zG £2, let 5^ G Uh be the smooth 5-function which satisfies. 


{^hAh) = Vh{z), zG Q.VvhGUh- 


For any appointed direction L, define the direction derivative 

dzv{z)= lim (v(z-l-Az)-v(z))/| Az|. 

Az-^0,Az//L 

Then, there exists a regularized Green function of derivative type (x) G //q (X2) f] (f2) 

such that 

a{v,dfi^) = (5^5;^, v), Vv G Hl{a). 

Let G Uh be the finite element approximation of such that 

a{vh,d^G^ - dfil) = 0, Vv/, G t//,. 


Clearly, we have 


a{vh,dzGl) = a{vh,d^G^) = (5^5;^, v/,) = 5jV/,(z), Vv/, G t4. (4.9) 

Assume that partition 7), is quasi-uniform so that the inverse inequality holds on fi¬ 
nite element space Uh- Under this condition, the following boundness estimates were 
given in II19II21II 

\\d,Gl\\i<Ch-'\lnh\z, ||5,G^J|i 4 <C|ln/i|, (4.10) 
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where C is a positive constant independent of z G Q.. 

Theorem 4.1 Let partition 7/j be quasi-uniform, and u and m/, be the solutions of 
problem (O and (EJ), respectively, u G Then we have 

Wn^u — Uh\\\,oo < C/z^|ln/!|||M||3_oo. (4.11) 

Proof From (14.91) . error equation (12.9b . Theorem 3.1 and (14.10b . we obtain 
d,{uh- nhu)(z) =a{uh-nhU,dfil) 

= a{uh - rifiU, dfil) - ah{uh - TlhU^nldfil) + atfuj, - TlhU^nldfil) 

= a{uh- TlhU, dfil) - ah{uh - TlhU^TIldfil) +ah{u- nhU,nj^dfil) 

< a{uh - TlhU^dfil) - ahiuh - nhU,nldfil) +Ch'^\\u\\i^^\\dfil\\\^\ 

< a(uh - TlhU^dfil) - ah{uh - TlhU^nldfil) +Ch^\\nh\\\u\\-i^o„. (4.12) 

Below we need to estimate (see Lemma 3.1) 

a{uh - rihu, dfil) - ah{uh - FlhU, 11*^3fil) 

= L I n-{A-A^)V{uh-nhu){d,Gl-nldfil)ds 

+ Y, [ n-A^v{uh-nhu){dfii-n;d,Gi)ds 

KeTi, 

+ ^ {-div{AV{uh-nhu)),dfii-nidfii)K 

KeTi, 

+{c{uii - Hhu ), dfil - nldfil)K 

= 51+52 + 53 + 54. (4.13) 

First, using (12.11b . (12.12b and (12.15b . Theorem 3.2 and (14.10b . we have 

5i < C ^ /i^|A|i,„||V(M,-Fr,M)||o.5i,||5,G;;-Fr;:a,GJJ|o.air 

KeT,, 

<ch^\viu,i-nhu)\\h^\dfii\u 

< Ch^\\uh\\dfilh <Ch^\\nh\^\\uh. 

Next, it follows from Lemma 4.2 and inverse inequality (12.12b that 

h^ 

Si = < Y, ni,u)xy\Q,oo,K\idzGl)xy\o,oo,K 

^ Y ^K\i‘^h — ^|,u)xy\Q^oc,K\\^zGl\\l,K■ 

KeT,, 















16 


T. Zhang, L.X. Tang 


From (I4.51 l and Lemma 4.1, we have 


\{uh — ni,u)xy\ = |(m/i — F//,m)i234| < V2\uh — ni,u\i^[c,h < y/^'^YWh — ni,u\i^K- 

Substituting this estimate into S 2 and using Theorem 3.2, we obtain 

S2<Ch^Uh-nhuh\\d,Gl\\i<Ch^\lnh\^\\uh. 

Similarly, by Lemma 4.1 and Theorem 3.2, we have (noting that (u^ — nhu)xx = 
{uh ^h^)yy ~ 0 ) 

Ken 

< C ^ \\u,-n,u\U^K\\d,Gl-n*,d,Gl)\\o,K 
Ken 

<Ch^\\uU\dfilh<Ch^\\nh\^\\uh. 

For ^4, using Theorem 3.2, (I2.141 i and (14.1 011 to obtain 

54<C ^ \A\i^^\\uh-n,M\i.K\\dzGl-md,Gl\\o,K 

Ken 

< C/!^||m||3/!||5^G^J|i <C/!^|ln/!|5||M||3. 

Substituting estimates Si ^ S 4 into (14.131 1. we complete the proof by (14.121 i. □ 

Now we can give the main superconvergence result. 

Theorem 4.2 Let partition 7), be quasi-uniform, and u and m/, be the solutions of 
problem (O and ( 12. (SI) , respectively, u S VT^’°°(X2). Then we have 

max|(VM — Wuf){P)\ < C/r^|ln/z|||M||3 <x>, (4.14) 

PeS 

where Vw{P) denotes the average values of gradient VwfP) in elements containing 
point P and S is the optima stress point set composed of all interior mesh points, the 
midpoints of interior edges and the midpoints of elements. 

Proof From (14.Il l and Theorem 4.1, we immediately obtain the conclusion by using 
the triangle inequality. □ 
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5 Numerical example 

In this section, we will present some numerical results to illustrate the theoretical 
analysis. 

Let us consider problem (ll.51 l with the data: 

[e^^+y'^ + l e^+y \ 

A{x,y)=\ , c{x,y)=2+x + y. 

y e^+y e^y+x^ + lj 

We take Q = [0,1]^ and the exact solution u{x,y) = 2 sin(27rx) sin(37ry). 

In the numerical experiments, we first partition domain Q into a rectangular mesh 
7/,, then the refined meshes are obtained by connecting the midpoints of each edge 
of elements in 7},. Thus, a rectangular mesh sequence is generated with successively 
halving mesh size /i/2', i= 1,2 ,---. 

Denote by e/, = max|VM(P) — Vuf,{P)\ the computation error and the numerical 
convergence rate is computed by using the formula r = ln(e/,/e/,/2)/ln2. Table 1 
gives the error and convergence rate with successively halved mesh sizes. We see that 
the convergence rate is just about 0 {h^), as the theoretical prediction. 


Table 1 Convergence rate of gradient approximation on set S 


h 

1/4 

1/8 

1/16 

1/32 

1/64 

1/128 

eh 

1.212 

3.099e-l 

7.856e-2 

1.969e-2 

4.949e-3 

1.243e-3 

rate 

- 

1.9671 

1.9802 

1.9961 

1.9926 

1.9932 
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